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We analyze the transverse momentum dependent distribution and fragmentation functions 
in space- like and time-like hard processes involving at least two hadrons, in particular 1-particle 
inclusive leptoproduction, the Drell-Yan process and two-particle inclusive hadron production 
in electron-positron annihilation. As is well-known, transverse momentum dependence allows 
for the appearance of unsuppressed single spin azimuthal asymmetries, such as Sivers and 
CO , Collins asymmetries. Recently, Belitsky, Ji and Yuan obtained fully color gauge invariant 

expressions for the relevant matrix elements appearing in these asymmetries at leading order in 
an expansion in the inverse hard scale. We rederive these results and extend them to observables 
at the next order in this expansion. We observe that at leading order one retains a probability 
interpretation, contrary to a claim in the literature and show the direct relation between 
^ i the Sivers effect in single spin asymmetries and the Qiu-Sterman mechanism. We also study 

fragmentation functions, where the process dependent gauge link structure of the correlators 
is not the only source of T-odd observables and discuss the implications for universality. 
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The study of polarization and transverse momentum dependent distribution functions was initiated by Ral- 
ston and Soper [1]. Their study of the Drell-Yan process was performed at tree level and did not address the 
color gauge (non-)invariance of the distribution functions. At leading order (tree level) no single spin asymme- 
tries were obtained in the Drell-Yan process (see also [2]). Sivers [3] proposed a specific non-trivial correlation 
involving polarization and transverse momentum, that would lead to unsuppressed single spin azimuthal asym- 
metries. For distribution functions such a correlation seemed to entail a violation of time reversal invariance. 
Collins [4] showed that this was not the case for similar correlations in the fragmentation process. Nevertheless, 
phenomenological studies of the consequences of the Sivers effect were performed [5,6]. Recently, Brodsky, 
Hwang and Schmidt [7] (BHS) demonstrated in an explicit model calculation that the Sivers asymmetry can in 
principle arise, after which Collins [8] demonstrated that it is the presence of a path-ordered exponential in the 
definition of transverse momentum dependent distribution functions that allows for the Sivers effect without a 
violation of time reversal invariance. 

This generated renewed interest in the proper gauge invariant definition of transverse momentum dependent 
correlators. The definitions of transverse momentum dependent parton densities of Ref. [9,10] did contain path- 
ordered exponentials (links) to ensure color gauge invariance, but these were not closed paths (each quark field 
has a straight link to infinity attached to it, but pieces at infinity are missing). If one includes such links by 
hand, it is no problem to consider also closed paths, but Efrcmov and Radyushkin [11] had demonstrated that 
the path of the link can be derived in transverse momentum integrated parton densities and this can also be 
done when the transverse momentum is not integrated over. In this way different processes can yield different 
paths [12], but no physical observable effects were expected from such links. However, until recently these 
derivations were incomplete, since the obtained paths were not closed. The missing piece would have to involve 
transverse gluon fields at lightcone infinity, which were thought not to affect physical matrix elements or at 
the very least lead to contributions suppressed compared to the leading order. Recently, the derivation of fully 
color gauge invariant matrix elements, with paths closed at light-cone infinity, was completed by Belitsky, Ji 
and Yuan [13,14]. They observed that transverse gluon fields do not always lead to suppression, contrary to 
common belief, formalizing the model results of BHS. The resulting fully color gauge invariant matrix elements 
strengthen the observation of Collins [8] that the presence of the link invalidated the earlier proof of the absence 
of the Sivers function due to time reversal invariance. 

With all these technical details clarified, the justification of the phenomenological studies of the Sivers (and 
similar) effects was provided. Next, however, the question of observable process-dependence arose. Collins [8] 
demonstrated that the Sivers asymmetry in (semi-inclusive) deep inelastic scattering (DIS) and the Drell-Yan 
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process must occur with opposite signs. This has been confirmed in the BHS model calculation [15], but still 
awaits experimental verification. It would be the first demonstration of an observable effect due to the presence 
of a path-ordered exponential in the hadron correlators and thereby would show the intrinsic non-locality of the 
operators occuring in these semi-inclusive processes. 

Given this process dependence it is relevant to study which color gauge invariant distribution and fragmen- 
tation functions appear in different processes. In this paper we study color gauge invariance of transverse 
momentum dependent distribution and fragmentation functions appearing in hard processes, in particular in 
semi-inclusive deep inelastic leptoproduction (SIDIS), the Drell-Yan process (DY) and e + e _ -annihilation. 

We will employ a field theoretical approach to these hard processes and follow the notation and derivation 
of Ref. [12], now taking into account the additional contributions uncovered by Belitsky, Ji and Yuan [14]. 
Our analysis is different at several points, but we confirm their results. In addition, we obtain new results 
for the first sub-leading order results in an expansion in inverse powers of the hard scale. For instance, we 
demonstrate for the first time a direct relation between the Sivers effect in single spin asymmetries and the 
Qiu-Sterman mechanism. Also, we study transverse momentum dependent fragmentation functions, where the 
process dependence of the gauge link structure of the correlators is not simply an overall sign. Rather one finds 
that two different (but universal) matrix elements enter in different combinations. 

In this paper 'leading' and 'sub-leading' always refer to the expansion in inverse powers of the hard scale. 
Perturbative QCD corrections beyond tree level (next-to-leading order in a s ) will need to be taken into account 
as well in further studies (see Refs. [10,4,16] for discussions of additional complications beyond tree level). In 
this paper we present a new way of isolating the leading and first sub-leading order parts of the cross sections 
in terms of correlators including the proper gauge links before evaluating them explicitly. These separate color 
gauge invariant expressions for each order have not been presented before. They facilitate the evaluation of 
asymmetries arising at a given order. 

Now we will outline the more technical steps to be followed in this paper. In the field theoretical approach, 
expressions for the structure functions of inclusive deep inelastic scattering (DIS) are obtained from diagrams 
as shown in Fig. 1 [17,18]. In these diagrams soft parts appear that represent matrix elements of the fields 
corresponding to the quark and gluon legs connecting the hard and soft parts of a diagram. The expressions for 
SIDIS structure functions are obtained from diagrams as shown in Fig. 2. Again soft parts represent specific 
matrix elements. In this paper we will only consider tree-level results, which means that if gluons appear, they 
are in essence legs of the soft parts. In other words, their (soft) couplings to the hard scattering part are included 
in the definition of the matrix elements. This approach requires a careful treatment of the diagrams involving 
quark-quark-gluon matrix elements such as those in Fig. lb or Figs 2b - 2f. This is important in order to arrive 
at color gauge invariant matrix elements that form the universal quantities, the distribution and fragmentation 
functions, appearing in cross sections. 

In section II we outline the diagrammatic approach in a number of steps, using two complementary lightconc 
directions n + and rt_, which in the presence of a hard scale (in DIS or SIDIS, the photon momentum) are fixed 
by the hadron momenta. The simplest (handbag) diagrams in Figs la and 2a only involve quark-quark matrix 
elements. In DIS the hadron momentum defines the lightcone direction n + and the nonlocality in the matrix 
elements is restricted along the lightcone direction n_ (for which n+ • n_ = 1). As is well-known, diagrams as in 
Fig. lb with any number of A + — A ■ n_ gluons yield the necessary gauge link connecting the two quark fields 
[11]. The nonlocal quark-quark operator combination with a gauge link can be expanded into a tower of local 
twist-two operators with different spins. Their matrix elements appear in the cross section as leading terms 
in an expansion in inverse powers of the hard scale. Diagrams with (transverse) A" gluons or with A~ gluons 
appear in matrix elements of higher twist operators, which appear in the cross section in terms suppressed by 




(a) (b) 
FIG. 1. Diagrams contributing in inclusive deep inelastic scattering 
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FIG. 2. Diagrams contributing in 1-particle inclusive deep inelastic scattering 



inverse powers of the hard scale. 

The situation in SIDIS (Figs. 2), discussed in section III, differs in a subtle way from that of DIS, because 
the nonlocality in the operator combinations is not restricted to the lightcone, but involves also transverse 
separations. The kinematics only constrain the nonlocality to the lightfront. In our analysis we first consider 
the A + gluon legs in diagrams as in Fig. 2b and 2f. These diagrams, as in DIS, will give rise to gauge links, 
but in this case connecting a quark field along the n_ direction to ±oo (cf. [10,12]), where the sign depends on 
the type of process. By including diagrams of the type in Fig. 2c as well, one can absorb all A + gluons into 
the lower blob. Diagrams like Figs. 2d and 2e allow one to absorb all A~~ gluons into the upper blob, resulting 
again in gauge links running along n + to infinity. Effectively one then considers only Fig. 2a, but now with a 
$ and A that contain the gauge links. 

Diagrams with transverse ^4"-legs, lead to quark-quark-gluon matrix elements, which will turn out to be 
suppressed, except for the boundary terms at lightcone infinity, recently discussed by Belitsky et al. [14]. We 
outline an alternative for this procedure and show that the latter appear when one expresses these fields in the 
appropriate field strength tensor G +a . The boundary terms that arise in this way combine into the transverse 
piece that completes the gauge link connecting the two quark fields (running via lightcone infinity). Upon 
integration over transverse momenta the result reduces to the correct gauge invariant operator of ordinary 
inclusive DIS. 

In sections IV- VI a comparison is made between different processes involving at least two hadrons, in particular 
between 1-particle inclusive leptoproduction, the Drell-Yan process and two-particle inclusive hadron production 
in electron-positron annihilation. For instance, in Drell-Yan the links run in opposite directions along the n_ 
direction compared to SIDIS (as noticed in [12,8]). Because the two situations can be connected via a time 
reversal operation, one can define T-even and T-odd functions that appear in the parametrization of the color 
gauge invariant matrix elements. For these functions factorization in principle should hold, although they 
appear with different signs in SIDIS and DY [8]. The T-odd functions appear in single spin asymmetries in 
these processes [3,4,19,6] or they appear in pairs in unpolarized azimuthal asymmetries [20-22]. In section 
VII we study the time reversal properties of distribution and fragmentation functions and present explicit 
parametrizations . 
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II. HADRON TENSOR AND CORRELATORS IN SIDIS 



The hadron tensor for 1-particle inclusive leptoproduction is given by 

d 3 P x 

' X 

with being the product of current expectation values 

H^ H \P x ;P,S;P h ,S h ) = (P, S\J„(0)\P x ; P h S h )(P x ; P h S h \J„(0)\P, S). 



1 f P 

2MWl H Hq;P,S;P h ,S h ) = j^ J (2 ^ )32 po (2^ 4 (g + P Px P h )H l »{P x -PS; P h S h ). (1) 



(2) 



Due to the fact that the summation and integration over final states is not complete, prohibiting the formal 
use of the operator product expansion, we proceed along the lines of the diagrammatic approach of Refs. [1,17], 
based on nonlocal operators. The quark and gluon lines connected to the soft parts represent matrix elements 
of (nonlocal) quark and gluon operators. 

The hadron tensor is calculated for current fragmentation in deep inelastic scattering. In that case the 
exchanged momentum — q 2 = Q 2 is large and one has for the target momentum P and the produced hadron 
momentum P h the conditions that P ■ q, Ph- q and P ■ P h are large, of 0(Q 2 ). One is able to make a systematic 
expansion in orders of 1/Q, of which we will only consider the first two terms, (1/Q)° and (1/Q) 1 . In this 
situation one uses the scaling variables 
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where the approximate sign indicates equalities up to 1/Q 2 (mass) corrections. It is convenient to introduce 
lightlike vectors n + and n_ satisfying n + ■ n_ = 1 along the hadron momenta writing 



P^ = 
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(5) 
(6) 

(7) 



with q 2 = — Q 2 and Q 2 — Q 2 + Q 2 T . These equations define the lightcone coordinates a ± = a ■ n T and 
the transverse projector gf^ v = g^ v — vS+n} . In our treatment of the 1-particlc inclusive process we will 
consider Q 2 <C Q 2 , hence Q 2 w Q 2 , while £ w x B and Q w Zh up to mass corrections of order 1/Q 2 . The 
vector q£ ss + xP^ — P£ / ' z determines the off-collinearity in the process. In principle, mass corrections 
can straightforwardly be incorporated. Important to note is that the lightlikc directions n± = n±(P, Ph) are 
determined by the hadron momenta P and Ph- 

From the diagrammatic expansion (see Fig. 2a-e) one obtains up to 0(g), 
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FIG. 3. Soft parts representing the quark-quark and quark-quark-gluon matrix elements used in Eqs. 9-12 



where 

*y(p;P,S) = J 0^ e^(P,S\^(0)MO\P,S), (9) 

Ay(fc; Ph, Sh) = E/ (^I elfc ' 5 < 0|^(O|A,A- X P h ,X\^ 3 (0)\0 >, (10) 

*5 i >'P-Pi; P " S ) = / ^^^e** 1 ^- (^51^(0)^(^(01^5), (11) 



(2tt) 4 (2tt) 4 

A% ij (k,k-k 1 ;P h ,S h ) = Y / J ^-^e ik ^e ik ^-^(0\^(09A a (ri)\P h ,X>< P h ,X|^.(0)|0). (12) 

illustrated in Fig. 3. In the above expression we have omitted the contributions with the opposite direction on 
the fermion line. It adds to the result in Eq. 8 terms with q «-> — q and fx <-> v. In cross sections it will always 
lead to extending a sum over contributions from quarks to the sum over quarks and antiquarks. 

The aim of the calculation is an expansion in powers of 1/Q. For this a number of considerations are 
important. First, the matrix elements represented by blobs in the diagrammatic expansion should vanish fast 
enough when any of the products of momenta involved becomes large, e.g. the virtualities of the quarks or 
gluons. To be precise, in Fig. 3 the products p 2 ~ p\ ~ p ■ pi ~ p ■ P ~ p\ ■ P ~ P 2 = M 2 <C Q 2 ■ With 
the choice of parametrization in Eqs. 5-7, this implies that for the momenta in Figs 3a and 3c one has for the 
plus-components p + , p^, P + ~ Q, while the minus-components p~ , p± 7 P~ ~ 1/Q. For the fragmentation parts 
(Figs 3b and 3d) one has minus-components k~ , fc-f, P// ~ Q, while for the plus-components k + , kf, Pj[ ~ l/Q. 
The transverse momenta are of 0{M). Introducing momentum fractions x = p + / P + and z — Pf~ /k~ , writing 

p^ = p- + x P+ n% + p£, (13) 

k" = ^-n^+k + n'-l + k^, (14) 

one finds, when neglecting 0(1/Q 2 ) contributions, that 8 A {p + q — k) — ► 5(x — x B ) 5(z — zy/) S 2 (p T + q T — k T ), 
thus identifying the scaling variables and momentum fractions, x — x B and z = Zy\. 

Thus, in a calculation up to 0{1/Q 2 ), the integration over the minus-components of momenta in the matrix 
elements $ and $a can be performed, restricting them to the lightfront, 



*ii(?,Pr) = J dp- $ntp;P,S) = J dg (2 ^ T jP< < P,S\^(0)MO\P,S> 



(15) 



while &a(p + 7 Pt , pf , Pit) = f dp dp 1 $2 {pi P ~ Pi i P 5) involves two integrations over the minus-components 
of the parton momenta. We will occasionally also use the variable x\ defined via p\ = x\ P + . The integrations 
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over minus-components are sufficient to render the time-ordering in these matrix eiements superfluous, which 
can be proven completely analogous to the proof for the matrix elements in which also the integration over 
transverse momenta is performed, in that case restricting them to the lightcone [23,24]. In the matrix elements 
of the types A and A^ the integrations over the plus-components of the quark and gluon momenta can be 
performed, leading to lightfront correlation functions A(z, k T ), 



A y -(z,fc T ) = J dk+ A l3 (k-P h ) = E/ ^ < °\MO\Ph,X >< Pfc,X|^.(0)|0 > 



(16) 



and A°[(fc-,fc T ,fcf ,k lT ) = J dk+ dk+ A^(fc, k - h; P h ). 

Matrix elements like <E>(x,p T ) have a particular Dirac structure, Lorentz structure, and canonical dimension, 
which must be visible in the parametrization of $ through the dependence on non-integrated parton momenta 
and the hadron momentum and spin vectors. One deduces immediately the Dirac structure J dp~ $ ~ fi + = 7 ~ 
giving a leading contribution and the Dirac structure involving the unit matrix in Dirac space requiring in 
addition a factor P~ in the parametrization, which will lead in the calculation of Eq. 8 to a suppression factor 
1/Q. The leading structure of J dp~ dp^[ $^ matrix elements involving two integrations over minus components 
gives for J dp~ dp^[ <I>j[ ~ 7 ~ , similar to J dp~ <£>, but for a transverse gluon one gets J dp~ dp^ ~ P 
leading to a 1/Q suppression in the calculation of Eq. 8 (apart from the subtlety with the boundary terms, 
where the role of P~ is taken over by S(pf), to be elaborated upon below). Of course in a parametrization of the 
latter matrix element also the transverse index must appear, e.g. a non-integrated parton transverse momentum 
p T or the spin vector S T in case of a transversely polarized hadron, but these are not relevant for an expansion 
in powers of 1/Q. The matrix element J dp~ dp~[ $^ will always appear suppressed by at least (P~) 2 — ► 1/Q 2 - 

For the fragmentation parts one has after integration over plus-components J dk + A ~ J dk + dkf A^ ~ 7+ , 
while f dk + dk^/ AJ T <~ P^ 7+, becoming suppressed by 1/Q (again apart from the subtlety with boundary 
terms). The matrix element J dk + dkf A^ will always appear suppressed by at least (P^) 2 — > 1/Q 2 - The 
explicit parametrizations for the matrix elements in terms of distribution and fragmentation functions have 
been extensively discussed in many papers [25,2,19,6,20] and will be summarized in section VII. 

In order to find the leading contributions, we need in the calculation of Eq. 8 not only the first term (diagram 
in Fig. 2a), but also the terms involving $^ and A^ (Figs. 2b-e) and even multiple-gluon matrix elements of the 
form (Fig. 2f), etc. Such a resumming of multiple-gluon matrix elements can be easily performed in DIS, 
where the integration over transverse momenta of partons can always be performed in addition to the minus- 
integration. The resummation leads to a modified first term in Eq. 8 with in the &(x) — J dp" d 2 p T $(p; P, S) 
matrix element the inclusion of a gauge link 



tfce = p <■*!>( -hi T d C^ + (C) 



EHfff -1 U,'... / d( N A+(( N ) 



(17) 



connecting the quark fields, rendering the object color-gauge invariant [11]. We will discuss the full procedure 
to obtain a color-gauge invariant object in 1-particle inclusive leptoproduction in the next section, following in 
part Ref. [12] and recent work by Belitsky, Ji and Yuan [14]. 

III. COLOR GAUGE INVARIANCE IN SIDIS 

In this section we will discuss the resummation of contributions in SIDIS coming from diagrams in Figs 2b-f. 

At orders (1/Q)° (leading) and (1/Q) 1 (first sub-leading) the integrations overp - , p^ , and k + in the corre- 
sponding soft parts can be performed. The result of the first term of four quark-quark-gluon contributions in 
Eq. 8 is 

[term 1] = - J d 4 p d 4 k 5\p + q - k)^J d 4 Pl Tr ( 7a - - 7y $S(P,P - Pi) 7/ ,A(fc)) j 

= -J d 2 p T d 2 k T 8 2 (p T + q T - k T ) dp+ d 2 PlT j ^ e'^e i "'-("-fl 
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TABLE I. The color gauge covariant objects into which the gluon fields in SIDIS are combined depending on the Dirac 
structure of specific terms in the hard quark propagator to which the gluon couples 
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G +a 




G +a 




U T 



x (p, s\m^(z, k T ) la {k J t p J 1 _ + T ^ +te i, 9 A a (v)m\p, s) 



(18) 

£+=r,+ =0 



where is made explicit (Eq. 11) and the minus- and plus-integrations are performed. In the expression 
after the second equal sign, it is understood that in the integrand p~ = = k + = while p + = x P + and 
kr = P-/z. 

Next we will split off from the quark propagator those parts that are relevant at leading and first sub-leading 
order in l/Q. These parts depend on whether the index a of the gluon field is plus, transverse or minus. For 
the l/Q order, we will restrict ourselves to obtaining a color gauge invariant expression for the hadron tensor 
integrated over the transverse momentum q T of the photon. For this result one first needs to consider the 
leading order term unintegrated over q T . The end results for the hadron tensor in several cases are summarized 
in the next section. 

One finds for the quark propagator explicitly (with k~ q~ = Q/\/2) 

fi-fa+m _ {ft + rn) -fi+pf ~i>\T ^ 



(k - pi) 2 - to 2 + ie - p + QV2 + (k T - p lT ) 2 - m 2 + ie 

Obvious contributions at leading order are, the k~fi- term of the quark propagator, which in combination 
with A + gluons leads to the link operator in the rf~ direction. Less obvious are the contributions from fields 
that are independent of rf~ which, as can be seen from Eq. 18, lead to a delta-function S(pi). In that case other 
leading contributions appear. In particular, a contribution coming from the last (transverse) term will lead to 
(leading) link contributions in the transverse direction. 

Contributions at order l/Q are coming from the fi + term of the quark propagator in combination with the 
transverse gluons, and the transverse part of the quark propagator, $i T , in combination with the A + gluons. 
These contributions can be combined into a color gauge invariant matrix element containing the field strength 
tensor. A summary is given in Table I. 

The leading contribution in Eq. 18 comes from <I>^. We use that 7~($ + m) = 2k~ — (]t — m)7~, the fact 
that A(fc) (^ — m) ~ <~ l/Q (QCD equations of motion) to obtain 

A(fc)7 -f A + (n) W -A(fc)7 y-r tA + (t}) 

(k - Pl ) 2 - m 2 + ie K " W ' 2k-(p+-ie) " 

« _A(jfe) A+{V) + A(k) Ijtl^M + A(fc) (l/- m ) 1 -A+(r 1 ) 
p\ — ie 2k~ (pi — ie) 2k~ (pf — ie) 

--A( fc) ^ + A (fc )^^M (20 ) 
pt-ie QV2 (pt - ie) 

with omitted parts being of Oil/Q 2 ). The first term inserted in Eq. 18 gives a leading contribution, 
[term 1.1]= | d 2 Pr d 2 fc T 5 2 (p T + q T - k T ) dp+ d 2 p lT J d ^.^ e ^ e ^<v-0 



x (p,s\m^(z,k T ) lv ?^-ip(o\p,s) 

Pi - *e 



J dV d 2 fc T 5 2 (p T + q T - k T ) J ^^e*^ 

x <P,S\^(0)^A(z,k T h v (-ig) f dry- A+(r,)^)\P, S > 



(21) 



£+=»)+=0: Vt=£t 
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This is precisely the O(g) term in the expansion of ^ multiplying The result of the diagram in Fig. 2f 

with two ^4+-gluons gives the 0(g 2 ) term, etc. From the second term in Eq. 18 (diagram in Fig. 2c) one obtains 
the 0(g) term in the expansion of ^ following ip(0). 

The A^-gluons in the other diagrams in Fig. 2d and e and corresponding higher orders can all be absorbed 
into link operators in modified soft parts of the form 



dT d 2 £ r 
(2tt) 3 

d£+ d 2 e T 



> 



(22) 
(23) 



where [/^^ indicates a link along the lightcone plus-direction running from — oo to £ + . These quantities, 
however, are not color gauge invariant, although we note that upon integration over p T and k T one obtains 
color gauge invariant lightcone correlators <&(x) and A(z), in which the two links merge into one connecting 
the lightlike separated points and £. These are e.g. important in q T -integrated SIDIS cross sections at leading 
order. For the transverse momentum dependent functions, however, we are still missing a transverse piece 
that leads to color gauge invariant definitions. It has to come from transverse gluons, which are next to be 
investigated. 

Since the dominant part of j dk + A is proportional to 7+, one finds (naively) for the transverse gluons in 
term 1 (Eq. 18), 



^ 1 + m ^)«A( %tt -^-i ?W 



(k — pi) 2 — m 2 + ie 



-A(fc) 
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QV2 



(24) 



The (remaining) second term in Eq. 20 and the result of Eq. 24 give as 0(1/ Q) contribution in term 1, 

d£~ d 2 g T d?T dV ip< 



[term 1.2] = / d p T d k T S (p T + q T — k T ) dpj~ d p\ T 
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J d 2 p T d 2 k T S 2 (p T + q T — k T ) d 2 pi T J 
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^= <P,5|^(0) 7M A(z,fc T ) 7 - 7Q IaZ(£)-J' dn-d?A+(ri-) \ 7 ,V(0I^5> 



QV2 



(25) 



Including in addition all diagrams with longitudinal A + gluon fields, all colored fields become linked along the 
minus-direction, with the same link directions for $ and $,4- Using the relation between G +a and A^, outlined 
in the Appendix A including all minus links, we find (suppressing the links U~) 



K(0- I drT (a?A+(rT))= / d v - G +a ( V ) + A«(oo-), 
with the points r]~ = (t]~ , , £ T ) and oo~ = (00, £+, £ T ). The part of term 1.2 containing G +a is 
[term 1.2a] = J d 2 p T d 2 k T 5 2 (p T + q T — k T ) J ^ 



(26) 



Jp-Z 
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<P, S\m 7Al A^,fc T ) 7 - 7Q | 4 dr G +a (f 1 ) lu m\P,S 



QV2 



> 



(27) 
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Upon integrating the hadron tensor over transverse momenta q T , the above convolution factorizes and produces 
a color-gauge invariant 0(1/Q) term, 

J d 2 q T [term 1.2a] = Tr ( 2_| lv J°° dp + _l_.^(p+,p+-p+)^A(z^j , (28) 
where (including link operators) 



(29) 



also denoted &q(x,x — Xi) with X\ = p^ /P + , and with LC denoting {£+ = i] + = S, T = Vt = 0}. 

We are left with a boundary term containing A T (oo~), which needs special care. The argument of the 
transverse field in the boundary term is fixed by the link direction in U~. The consequence is that the rf~ 
dependence disappears. We note that the integration over r\~ thus can simply be performed, showing that one 
deals with a matrix element that is proportional to S(pi) in momentum space [26], 

*(tf)*A(oo)«(P.P-Pi) = / c i ^e i ^-("-«<P,5|^.(0)^(oo,f, + ,f h .)^(0|P,S). (30) 

Because <5(p^~) ~ 1/Q one finds that J dp" dp\ $2(oo) ~ 7 + > 1 - c - ^ ^ s not suppressed. This means we have to 
revisit the approximations made to the fermion propagator for the boundary term. Going back to the starting 
point in Eq. 18 we obtain for the boundary contribution after integration over r)~ , 

[term 1.2b] = -J d 2 p T d 2 k T S 2 (p T + q T k T ) dp+ d 2 PlT J ^~^ T e ^V^^-^> S{p +) 



x < P, S\ ^(0) 7ai A(z, k T ) la - * . 7^(00, v + ,Vt) m\P, S > 



£+=,,+ =0 

(31) 



After a substitution for A T , 



f r]T 

g AZ( V ) = i8%(-ig) d( T -A T (ri-,V + ,( T ), (32) 

J COT 

we do a partial integration. In the matrix clement, we then encounter the following part, which we need to 
consider in the soft gluon limit (pf = 0) in which the denominator of the quark propagator can no longer be 
approximated as in Eq. 20 or 24. Realizing that p\ <~ Q and p~[ <~ 1/Q we obtain 

-At*)(-» + A + ») (t _^*_ + J +fc < W) 

= -A(*)«(p+). (33) 
The result after integration over i] T , p\ and pi T is a term 

q2, 



[term 1.2b] = ^ d 2 p T d 2 fc T <5 2 (p T + g T - fc T ) ^ ^^le^ 

x <P, S\^(0)^A(z,k T ) lu (-ig) f T d( T -A T (^,OXr)ip(0\P,S> 



(34) 



which gives precisely the first term of the transverse link that is needed to modify Eq. 22 into a fully color gauge 
invariant matrix clement. Note that we did not assume a specific pure gauge expression for the At field at 
(~ = oo. Furthermore we did not neglect the quark masses in the quark propagator. For N transverse gluons 
we have to work out the boundary terms with more transverse gluons, for which we need the following relation 
that also holds for nonabelian fields, 
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i- 
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(a) 




GO 



(c) 

FIG. 4. Link structure for $' + '(a;,p T ) (a and b) and for &(x) (c) after integration over transverse momenta. 



( Hs) M f" dCiT • MCi) ■■■ T dCAf r • a t (Cm) ) a A%{v) 

= ia% ((-i9) M+1 /"dCiT • MCi) ... T cICm+it • ^t(Cm+i)) 



(35) 



This indeed produces nicely the nested integrations needed in the path-ordered exponential and we find that 
the terms 1.1 and 1.2b, are taken care of by using in the O(g ) part of Eq. 8 the color gauge invariant matrix 
element [14] 



dr d 2 e 



P,S> 



(36) 



with 



V, 



V exp y dCr • A T (() 

jt r d ci- ■ mci) ■ ■ ■ r dCAr r • ^ 



C+=C+=Q+,C-=?-=a- 



■(Cat) 



(37) 



The (7jq t ^^j-link of course comes from diagrams with transverse gluons like in Fig. 2c and higher orders. The 
link structure for the soft part describing the distribution of quarks in a hadron probed by a spacelike photon 
is illustrated in Fig. 4. The direction of the link, running to +oo along the minus direction is indicated via the 
superscript [+] in Eq. 36. In other processes one will find that the link can also run in the opposite direction to 
—oo along the minus direction. This will be indicated with a superscript [— ]. 

Eq. 36 is an important expression, since in its full generality it allows for certain distribution functions, 
usually referred to as T-odd functions, that would be absent in case one ignores the gauge links (cf. e.g. [3,4,6]). 
Without the transverse gauge links, it may therefore seem that a choice of A + = gauge would demonstrate 
the absence of such T-odd functions. In the derivation of Eq. 36 no gauge was assumed (one can actually arrive 
at this result by first considering the A~ — gauge as done in Ref. [12]), hence it should not be viewed as one 
out of many ways to "gauge-invariantize" the matrix element [27]. Up to the order we consider here, the result 
is derived rather than assumed. 

We note that in leading results the color-gauge invariant object $[ + '(x,p T ) contracted with 7" 1 " still is a semi- 
positive definite matrix in Dirac space, which is the basis for deriving positivity conditions, such as the Soffer 
bound [28] and many more [29]. Hence, we disagree with the statement "Structure functions are not parton 
probabilities" by Brodsky et al. [30] . To be precise, the distribution functions containing the transverse link are 
still probability densities. 

As mentioned already, when one considers the g T -integrated results one obtains the lightcone quark-quark 
correlations with the link in Eq. 17. The transverse link does not affect that result and one has & + \x) = 
i>[ _ l(a;) = (see Fig. 4c). If one looks at azimuthal asymmetries or weighted cross sections one needs to 

consider matrix elements weighted with transverse momentum. In those cases one explicitly needs to take into 

account the transverse part of the link. We define transverse moments $g a (x) , 



(38) 



10 



that in a straightforward way can be related to color gauge invariant quark-quark-gluon matrix elements 
and $2>: the latter involving the covariant derivative, 



(P,S\^,(0)U m iD^)\P,S) 



LC 



-<P,S|^-(0)^- 



±00] 



±00 



LC 



or 



where 



/oo 
-00 Pi T »e 



*£«(*) = / $s«(p + ,p + -pr) 



^1 

2tt 



e+ = o 



(39) 



(40) 



(41) 



LC 



The important observation we want to make here is that the difference between correlation functions with links 
running to ±00, respectively, is related to a quark-gluon correlator, 



(42) 



the latter being given the name gluonic-pole matrix element since it corresponds to the soft-gluon point = 0. 
Its consequences have been studied for several processes [31-35,26] and it is viewed as one of the possible 
mechanism to generate single spin asymmetries. We will comment on this further below, but already mention 
that the above relation between <f>^ a (x) and $q(x, x) implies that the Sivers effect [3] is directly related to 
the Qiu-Stcrman mechanism (the gluonic-pole matrix clement), i.e. if one is nonzero, then the other also is. We 
will make this relation more specific below. 
We further define 

n(x)^l(^ ]a (x)+^> [ d ]a (x)), (43) 

$2(a;) = PV / d Xl — x - xx), (44) 
J xi 

where we use &a(x) to distinguish the function from the non-gauge-invariant $J(x). These definitions imply 

$k ±la (s) = $g(s)±7r (45) 



n{x) = $ a D {x)-3> a A {x). 



(46) 



The relations in Eqs 44 - 46 are relations connecting color gauge-invariant quantities. We will return to the 
above functions and their properties in section VII. 

We end this section by giving the remaining contributions in the g T -integrated hadron tensor at order 1/Q. 
We give a systematic summary of the SIDIS hadron tensor for several cases in the next section. We can use 
Eq. 40 to rewrite the twist-3 contribution in Eq. 28 obtained after ^-integration into the form 



J d\ T [term 1.2a] = Tr ( 2_| 7 „ [$« {x ) - d>W»] ^ A(*)) . 



(47) 



The transverse gluons in the second term of Eq. 18 (diagram in Fig. 2c) produces besides the transverse link 
£/jo t00t ], already absorbed into Eq. 36, also a twist-three piece. For this one needs matrix elements with 
interchanged arguments such as $a(p — Pi,p)- The resulting twist-3 term after integration over transverse 
momenta is 

/ d\ T [term 2.2a] = Tr ( lv f°° dp+ -p— <&g(p+ - pt,P + ) 7m 7^7= A(z) 
J V J -00 Pi + *e QV2 



Tr 7, 



Pi + it 
7o*S + (»)7o -4 +la W 



7a 7 a / ' 



QV2 



(48) 
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where we have used that <$ffl(x) — 70 &^ a ^(x) 70. The third and fourth 0(g) terms in Eq. 18 (diagrams in 
Fig. 2d and e) need the (direct and conjugate) quark propagators (with p + w — q + = Q/y/2) 

•p — fci + m (j/> + m) — fi- ki — fci T 



[p - fci) 2 - to 2 ± ie QV2 + (p T - k lT ) 2 -m 2 ± 



(49) 



The calculation yields the links « and running along the plus-direction in the matrix elements 

A(z,k T ) in Eq. 23 for the fragmentation part. Including the transverse gluons we get the fully color-gauge 
invariant matrix element indicated as the spacelike fragmentation A\~\z,k T ) with the link as indicated in 
Fig. 5. Furthermore one obtains twist-3 contributions containing Ag(fc — ki,k) and Ag,(fc, k — k\), which after 
integration over q T yield 



d q T [term 3.2a] = Tr I 7^ 



dk; 



QV2 



Tr 7/1 



7,^$(x)l, 



d 2 q T [term 4.2a] = Tr 



( 7 + 7a 

\Qy/2 



7m 



dfcf 



A£(*)-A Hq (*) IvHx) 



QV2 



A^(k~,k- -k^)^(x) 



(50) 



(51) 



In deriving these results we stepped over some subtleties involving the inclusion of diagrams with crossed gluon 
lines and the use of the full QCD equations of motion, but for this we refer to Ref. [12]. 

IV. SIDIS AND DIS CROSS SECTIONS 

The basic expression for W^ v (q; P, S; Ph, Sh) contains a convolution of the transverse momentum dependent 
functions. We have seen that upon integration over q T , 

J d 2 q T d 2 p T d 2 k T S 2 (p T + q T — k T ) . . . = J d 2 p T d 2 k T . . . , (52) 

the integral can be deconvoluted. This is also true for azimuthal asymmetries constructed by weighting with 

9?, 



J d 2 q T q% d 2 p T d 2 k T S 2 {p T + q T -k T )... = J d 2 p T d 2 k T (jfej - p°) . 



(53) 



If one calculates the 0(1/ Q) result, one has to be careful, however. One cannot simply perform the integration 
over q T in the hadron tensor, since the lepton tensor involves q. To proceed, one starts with a (Cartesian) set 



(a) spacelike distribution <E>[ H 









1 1 > 1 1 * < 




"1 




' i- 



(b) timclikc distribution & 











■ ' 


4 + 



(c) timelike fragmentation AM (d) spacelike fragmentation A^ 

FIG. 5. Link structure for $(x,p T ) (a and b) and A(z,k T ) (c and d). 
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of vectors, starting with q defining in SfDIS a spacelike direction while the other external vectors are used to 
define orthogonal directions. In particular, it is convenient to start with q and P = P — (P ■ q/q 2 )q, 

= ~ = -r, (54) 

* = 41 « + (55) 

with z 2 = — 1 and i 2 = 1. These two vectors can be used both for inclusive and semi- inclusive leptoproduction. 
The lepton tensor can be expressed in the vectors t and z and a perpendicular vector (equal for the lepton in 
initial and final state) which defines the azimuthal angle of the lepton scattering plane. From the Cartesian 
vectors two new lightlike vectors n'± = (t ± z)j\p2 = n'±(P,q) can be constructed as well as a perpendicular 

tensor g^ = g^ v + q^q v — tH v = g 1 *" — n^n'^K Since the lightlikc directions n'± are determined by P 
and q, instead of P and P h , the momentum of the produced hadron P h will have in general a nonvanishing 
perpendicular component enabling us to define a vector X^ = —P^j_/zh, which defines the azimuthal angle of 
the hadron production plane. 

It is straightforward to see that up to 0(1/Q 2 ) corrections, the previously defined set {n +1 ra_, q T } is related 
to the set {t, z, X} or {n' + , n'_, X} via 

w n'»-2^, (56) 
n% « n%, (57) 
<& « XV-yfi^n'*, (58) 

and X 2 w — Q^.. We note that the leptonic tensor is independent of X. The 1/Q term appearing on the 
righthandsidc of Eq. 58 is irrelevant in our calculations. Hence for experimental cross sections up to that 
order we can use for integrated and weighted cross sections the replacements J d 2 X . . . — » J d 2 q T . . . and 
/ d 2 X X a . . . — > J d 2 q T q™ .... We now consider separately integrated SIDIS, DIS and azimuthal asymmetries 
in SIDIS. We will also consider a few special cases. If one measures in the final state the jet-direction, this can 
be considered as a measurement of p T , i.e. q T = —p T . This is referred to as JET SIDIS. 

A. Integrated SIDIS cross section at leading order 

We have seen that inclusion of appropriate quark-gluon matrix elements make the O(g ) result in Eq. 8 
color-gauge invariant, 



2MW$(q;P,S;Ph,S h ) = J d 2 p T d 2 k T S 2 (p T + q T - k T )Tr ($M(x,p T )^A^(z,k T )^ . 



(59) 



At leading order the integration over transverse momenta of the produced hadrons simply can be performed 
and one obtains the basic result 



J d 2 X 2MW$(q;P,S;P h ,S h ) = Tr(*(x) ^ A(s) 7 „) 



n±—>-n 



± 



B. Azimuthal asymmetries in SIDIS at leading order 

We consider here cross sections obtained after integration over X and explicit weighting with X a . In practice 
this means measurement of the azimuthal angle of the produced hadron and compare it with other azimuthal 
angles, such as that of the lepton scattering plane, the (transverse) spin of the target hadron or the (transverse) 
spin of the produced hadron. For our purposes it implies calculation of J d 2 X I" 2MW M „, which at leading 
order gives 
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J d 2 X X a 2MWjS(q; P, S; P h , S h ) 

= Tr($(x) AH«(z) 7,) - Tr($[ +1 "(x) ^ A(z) 7 „) 



+ 



n± — >n. 

In these cross sections one finds for instance the Collins and Sivers effects [36,5,6,37]. 



QJ' 



(61) 



C. Integrated SIDIS cross section at 0(l/Q) 

As outlined one must be careful in integrating over transverse momenta. At 0(1/ Q) the differences between 
using n±{P,Ph) or n'±(P,q) matter. In particular the correlator A oc fi_ will lead to terms proportional to 
j^r/<3\/2- To find the ji- dependence in a Dirac space correlator, we use the projectors P± = 7 T 7 ± /2 = fi±fi T /2. 
We use that the leading term in A satisfies P_ A = A P + = P_ A P + to write 



A(z,k T ) = j (ft- fi+ A(z, k T ) fi + fi- 



A(z,k T ) 



1 / 



n± — * n± 



(jL T fi+ A(z, k T ) + A(z, k T ) fi + ^ , 



(62) 



and we obtain the 1/Q contribution coming from Eq. 59, 
J d 2 X 2MW$(q;P,S;Ph,S h ) = 0(1) result [Eq. 60] 

~07l/ d 2 q T j d 2 p T d 2 k T 6 2 (p T + q T - fc T ) j Tr(V+l (z, p T ) 7^ T ?4 + A^ (z , k T ) ^ 

+ Tr($[+1 (x,p T ) 7(U AH (z, k T ) fi+ i T 7 „) | 
= 0(1) result [Eq. 60] + Tr(^ ]a (x) 7, ^ A(z) 7,) - Tr(*(x) 7, ^ A [ d ]a (z) 7,) 
+ Tr(4 + 1>) 7, A(z) 0| 7,) - Tr(*(x) 7, A^ a (z) 0| 7,) • 



(63) 



Including these l/Q-terms and the four contributions from quark- gluon correlators with transverse gluons, one 
obtains the full integrated SIDIS cross section up to 0(1/Q), 



J d 2 X 2MW$ +1 \q;P,S;P h ,S h ) = 0(1) result [Eq. 60] 

+Tr ($M«(*) 7 , ^ A(z) 7,) - Tr($(x) 7, ^ A^fz) 7, 
- Tr (* 7, A(z) 0| 7,) - Tr($(x) 7, 7, 



+ ' 



+ Tr| y^ lv 



QV2 



~ *o ]a (x) 7m A(z) ) + Tr I 7, 7o^ T (*ho - <^ J » 



s Q t/ 



-J+l' 



7, ^ A(z) 



QV2 



+ Tr^ 7 p 
= 0(1) result [Eq. 60] 
+ Tr 



7oA«t( z)7o _AH«(z) 



7" 7 + ^ , rr, / 7 + 7a 



( h - 7, <TO 7m + Tr 7m A£(s) 



Tr 7m 4~ ]a (z) 7,*(*)) - Tr( 7 M ^ 4"^) 7, *(*)) + 0* - *)* , 



(64) 
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where the hermiticity properties of the various matrix elements have been used (see section VII). We note that 
the l/Q part in Eq. 63 is by itself not electromagnetically gauge invariant, but together with the parts arising 
from quark-gluon matrix elements, leading to the result in Eq. 64, it is gauge invariant. 

D. Integrated DIS cross section at 0(1/ Q) 

The hadron tensor for this case is obtained by integrating the SIDIS result over z and using for the fragmen- 
tation part the free (massless) quark — > quark result, A(z) = fi-S(l — z) and A^ a (z) — Ap(z) = 0. This gives 
the well-known result [17], 

2M W^(q; P, S) = Tr(<5(x) 7m 7 + 7 „) 

+ Tr ( 0| 7 „ n(x) 1, 7 + ) + Tr ( 0| 7 , 7v 7 + ) • (65) 

E. Azimuthal asymmetries in JET SIDIS at leading order 

Azimuthal asymmetries in JET DIS are obtained at measured q T = — Pj e t_L and with in addition still fixed 
P and q. We simply replace A(z, k T ) = fi- 5(1 — z) S 2 (k T ). Hence we start with 

2MW$(q;P,S;q T ) = Tr (*H (x, -q T )^7 + lJ) ■ (66) 

In this situation one will have to be careful with Sudakov effects [16], but the following weighted azimuthal 
asymmetry is free of these effects, 

J d 2 q T q% 2MWj?J(q; P, S; q T ) = -Tr (^ ]a (x) ltll + ^) . (67) 
This result is also a direct consequence of Eq. 61, taking for AH(z, k T ) the quark — ► quark limit. 

F. Other subleading cross sections 

Azimuthal asymmetries at 0(1/ Q), azimuthal asymmetries involving higher weighting than with one power 
of the momentum X a or SIDIS cross sections at 0(1/Q 2 ), require considerably more theoretical efforts than 
the one presented above. Moreover, it might be impossible to unambiguously disentangle hadrons within a 'jet' 
from hadrons in other 'jets', since the 'separation' of jets is merely an 0(Q 2 ) effect, meaning Pj c t i ■ fjet 2 <x Q 2 ■ 
This also implies that a factorization proof most likely cannot be given, just like the failure of factorization in 
unpolarized processes at l/Q 4 [38] . As is well-known, these difficulties do not appear for the inclusive DIS cross 
section involving just one (target) hadron which allows for a rigorous treatment at any order in powers of l/Q- 

V. THE DRELL-YAN CROSS SECTIONS 

For Drell-Yan, one has a similar treatment as for leptoproduction. The calculation involves now two soft 
distribution parts and annihilation of a quark- ant iquark pair into a gauge boson (we will only discuss the vector 
coupling here). The handbag diagram is given in Fig. 6a and an example of a diagram with an additional gluon 
in Fig. 6b. 

A full calculation at tree level including quark-gluon matrix elements as discussed for leptoproduction gives 
in this case 
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2MHV(g; P A , S a ;P b ,Sb) = / d 4 p d 4 fc 5\p + k-q){ Tr($(p) 7/1 $(fc) 7i ,) 



d pi Tr |^7 Q 
d 4 pi Tr I 7^ 



d 4 fcx Tr f 7v 
/ d 4 ^ Tr [ 7q 





— #6 — 


+ m 




(fc 


+ Pi) 2 - 


m 2 4 






-t-fr 


+ m 




(fc 


±Pi) 2 - 


m 2 — 


it 




i> + h- 


- m 




(p- 


ffci) 2 - 


TO 2 + 


it 




i> + h- 


f m 




(p 


+■ ^i) 2 - 


TO 2 — 


it 



la^>{ph^A( k - k l,k) 



+ . . 



(68) 



where $(p) and &a(p,P — Pi) are the same as in leptoproduction, but the role of A and A a is taken over by 

d 4 £ 



-ik-i 



(2tt)< 



(69) 



Aij (k,k- kx;PB,S B ) 



d 4 C dS . k 



(2tt) 4 (2tt) 



(Pb, ^(Off^M ^(0)|P B , 5 B ). (70) 



(note that this implies < &g(^, &t) = — fc a ^(^j k T )). 

The important difference between DY and SIDIS turns out to be the direction of the links. The result for the 
quark propagator in a quark-gluon diagram as in Fig. 2b, but then in the case of DY with a timelike outgoing 
photon (cf. Fig. 6a) yields a propagator 



— f, — f\ + m 



-m) - fi+pt - hi 



(k + Pl ) 2 ~ m 2 + it p+Qy/2 + (k T +pi T ) 2 -m 2 



it 



(71) 



(since k~ m q~ = Q/y/2). The difference with Eq. 19 is the sign with which pf appears in the denominator. 
For the j4 + -gluons this produces a link running along the minus direction to — oo, i.e. one finds transverse 
momentum dependent distribution functions &~^(xa,Pt), where x A ~ p + jP\ ~ q + /Pa - Also in the antiquark 
matrix element the link runs to — oo (along the plus direction with our choice of lightlike vectors), i.e. the A" 

gluons lead to the matrix element $ (xB,k T ), where x B w k~ / Pg ~ q~ / Pg . 

As in leptoproduction the tree-level calculation is most conveniently done with lightlike directions u±(Pa, Pb) 
defined via the hadron momenta. In order to perform the transverse integration at the level of the hadron 
tensor one introduces [39,40,2] a Cartesian set {t, z, X} starting with t = q/Q. A symmetric choice for z is 
z = (x A Pa — xb P B )/Q that defines the Collins-Soper frame. Using n' ± — (i ± z)/*/2, one has 



QV2' 



(72) 



4> 






q _ " 
























p il 




q 






(b 





FIG. 6. Quark-quark (a) and one of the quark-quark-gluon (b) correlators in tree-level diagrams for Drell-Yan scat- 
tering 
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while as in leptoproduction one finds that for the calculation of the cross section the orthogonal direction 
X differs from q T only by the irrelevant 0(1/Q) terms multiplying n'±, i.e. for our purposes X w q T rts 
q-X A P A - x B P B . 



A. Integrated DY cross section at leading order 



As indicated, the inclusion of appropriate quark-gluon matrix elements make the O(g ) result in Eq. 8 color- 
gauge invariant leading to 



W$(q;P A ,S A ;PB,S B ) = J d 2 p T d 2 k T S 2 (p T + k T - q T )Tr ($^(x A ,p T )^$ [ 1 (x B , fe r ) 7l/ ) . 



(74) 



This is analogous to the expression employed by Ralston and Soper [1], except that now the correlation functions 
are fully color gauge invariant. At leading order the integration over transverse momenta of the produced hadrons 
simply can be done and one obtains the basic result for the DY process, 



J d 2 X Wj?J(q;P A ,S A ;P B ,S B ) = Tr(<i>(x A ) lfl $(x B ) 7 „) 



+ ° 5'- 



(75) 



B. Azimuthal asymmetries in DY at leading order 



We consider here cross sections obtained after integration over X and explicit weighting with X a . In practice 
this means measurement of the azimuthal angle of the leptonic plane with respect to that of the hadron plane 
or relative to (transverse) spin azimuthal angles. For our purposes it implies calculation of / d 2 X X a W^v, 
which at leading order gives 



J d 2 X X a W$(q;P A ,S A ;P B ,S B ) 
= -Tr{^{x A ) 1 ^ [ d ]a {x B ) lv 



+ Tr 



(4 ]a {x A ) 1 ^{x B ) lv 



+ o 



Qj ' 



(76) 



C. Integrated DY cross section at 0(1/ Q) 



As outlined, above, one must be careful in the DY process in integrating over transverse momenta. In this case 
both correlators $ oc fi + and <f> oc fi- will lead to terms proportional to $ t /QV2- To find the fi- dependence 
in a Dirac space correlator, we again use the projectors P±. We now get 



$(x A ,p T ) = ^(fi + fi-§(x A ,p T )ifi-ifi^ 



<&(x A ,p T ) 



n± — > n 



~ T^J-f^ (ir i 1 - ®( x A,Pt) + ®{x A ,p T ) fl- i T ^j , (77) 



$(x B ,k T ) = ^fi-fi + ^(x B ,k T )fi+fi-^ 



$(x B ,k T ) 



- f l +^i x B,k T ) + $(x B ,k T )fi+ 4^. 



(78) 



Combining this 1/Q contribution coming from Eq. 74 with the parts from the quark-gluon diagrams one obtains 



17 




(a) (b) 



FIG. 7. Quark-quark (a) and one of the quark-quark-gluon (b) correlators in tree-level diagrams for back-to-back jet 
production in electron-positron annihilation 



J d 2 X wjZ +1) (q\PA,S A ;P B ,S B ) = 0(1) result [Eq. 75] 

- Tr ( 0| lv $«( Xa) 7m <p {xB ^J + Tr ^ W^(x B ) 7 **m) 

- \ Tr 7. *£~ 1q (sb) 7**(*a)) + 5 Tr( 7M ^ l^a*) 7, + (a* - *)* . (79) 
where the hermiticity properties of the various matrix elements have been used (see section VII). 



VI. BACK-TO-BACK JET PRODUCTION IN ELECTRON-POSITRON ANNIHILATION 

Also for 2-particle inclusive electron-positron annihilation we have a quite similar procedure. The calculation 
involves two soft fragmentation parts and the creation of a quark- antiquark pair. We will discuss only the case 
of creation from a (timelike) photon. The handbag diagram is given in Fig. 7a and an example of a diagram 
involving an additional gluon in Fig. 7b. 

The calculation of this tensor in a diagrammatic expansion proceeds as in the case of leptoproduction and 
gives 



W fl v(q\Pi,S 1 ;P2,S 2 ) = Jd 4 p d 4 fc 5\p + k-q)\ ■R(E(p) 7/ ,A(fr)7„) 



1 ; —pi / ~T~ PI H~ Tfl -r-Q' , . . / j \ 

d Pi Tr 7 Q -— ^ y— 7 l/ A j4 (p,p-pi)7 At A(fc) 

dV Tr ( 7^ — -7«A(fc)7vA^(p-pi,p) 
, ; , -r.. / - -f-h + m 



d fcl Tr A w Vrr 7aA(p) 7M A^(fc - fa,*) 

\ (p + fci)^ — + le 

where 

S y (p;P,5) =S / (04 e-^'«(0|^.(0)|P 2 ,X)(P 2 ,X|^(O|0), (81) 



x 
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A a Aij (p,p- Pl ;P,S) = ^J j^j^e 



d4 ^ d4 ^ e ~ip-t v-ipiin-t) 



(0\^JO)gA a (r 1 )\P 2 ,X}(P 2 ,X\Um)- (82) 



(note that Ag(z,p T ) = —p% A(z,p T )). 

It turns out that the direction of the links in the fragmentation functions changes in going from SIDIS to the 
annihilation process, just as the direction of the links in distribution functions changed in going from SIDIS to 
DY. The A~-gluons connected to hadron 1, produce a link running along the plus direction to +00, i.e. one finds 
transverse momentum dependent fragmentation functions AM (zi, k T ), where z\ s=s P^/k~ w P^/q~. Also in 
the antiquark matrix clement the link runs to +00 (along minus direction with our choice of lightlike vectors) , 

i.e. the inclusion of A + gluons lead to the matrix element A^ + \z 2 ,p T ) 1 where z 2 ~ P 2 /p + ~ P 2 /<? + ■ 

Again, to perform transverse integrations at the level of the hadron tensor one switches from n±(P\, P 2 ) 
directions to directions n'± = (t ± z)/\/2 with t = q/Q. In order to treat both 1-particle and 2-particle 
inclusive annihilation it is convenient [20] to fix z via one hadron momentum, for which we will choose P 2 , i.e. 
z = (q — 2 P 2 /z 2 )/Q. The relation between n± and n±, then is the same as in leptoproduction. The orthogonal 
direction determining the azimuthal angle of the hadron production plane, then is X = —P\±/z\ which for our 
purposes equals X m q T m q — P\jz\ — P 2 /z 2 . 



A. Integrated annihilation cross section at leading order 



The inclusion of appropriate quark-gluon matrix elements make the 0(5°) result in Eq. 80 color-gauge in- 
variant, 



WjS(q; Pi, Si;P 2) S 2 ) = Jd 2 p T d 2 k T S 2 (p T + k T - q T ) Tr (^A l+] (z 2 ,p T ) 1 ^+\z 1 , k T ) lv ) . 



(83) 



At leading order the integration over transverse momenta of the produced hadrons can be performed and one 
obtains the basic result 



J d 2 XWj°J(q;Pi,S i; P 2 ,S 2 ) = Tr(A(z 2 ) 7M A(z 1 ) 7 ,) 



+ O 



(84) 



By taking the free (massless) quark result A(z 2 ) = fi+ 5(1 — z 2 ) we obtain the 1-particle inclusive result (q T = 0), 



W$(q;Ph,S h ) = Tr(A(z h ) lul -l») + ° (^) ■ 



(85) 



B. Azimuthal asymmetries in the annihilation process at leading order 



We consider here cross sections obtained after integration over X and explicit weighting with X a , i.e. 
/ d 2 X X a W^u, which at leading order gives 



J A 2 X X a W$)(q; P\,S\:P 2 , S 2 ) 

= Tr(A(z 2 ) 1 , A [ + ]a ( Zl ) lv ) Tr(4 +lQ (^ 2 ) 1 , A{ Z1 ) lv ) 



O 



Q 



(86) 



By taking A(z 2 ) = fi + 5(1 - z 2 ) and Ag + '"(z 2 ) = 0, we get the weighted e + e — ► jet(Pj) + h(P h ) + X hadron 
tensor 



J d 2 X X a W${q;P h ,S h ;Pfi = Tr(A [ + ]a (z h ) 7,7-7,,) 



+ 



n±— >n± 



(87) 
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C. Integrated annihilation cross section at 0(1/Q) 



As before, one must be careful in integrating over transverse momenta at 0(1/Q). In particular the correlator 
A oc ?4_ will lead to terms proportional to ^ T /Q\/2. We obtain 



A(z\,k T ) = ^fi-fi + A(zi,k T )fi + fi- 



A(z u k T ] 



- ft+A(z 1 ,k T ) + A(z 1 ,k T )fi + 



(88) 



n± — > n ± 

which in analogy to leptoproduction leads to the full, integrated annihilation cross section up to 0(1/Q), 
J d 2 X W ( °+ 1 \q;P 1 ,S 1 ;P 2 ,S 2 ) = 0(1) result [Eq. 84] 

+ Tr ( 0| 7 „ A a D (z 2 ) 7 , A(z x )) - Tr 7 „ A£ ( Zl ) 7 ,A> 2 )) 

+ Tr (^7m A [ + ]a { Zl ) lv A{z 2 )^ -Tr( 7 , ^ A^^^ A> 2 )) + - „)* , (89) 

where the hermiticity properties of the various matrix elements have been used (see section VII). 

With the choice of i and z, the 1-particle inclusive result is found by taking A(z\) = fi_ 6(1 — z\), A^ a (z\) = 
A^(zi) = 0. After a (for ease of comparison) change of plus — > minus and taking the crossed (particle — > 
antiparticlc, [i ^> v) term we obtain 

W^,(g;P h ,S h ) = 0(1) result [Eq. 85] 

- Tr ( * AS( ^ >) - Tr ( > A ^'° ^) * • (90) 



VII. TIME-REVERSAL PROPERTIES 



In order to discuss the possible relations between the expressions for the various processes given in the previous 
three sections, we first discuss the time reversal properties of the correlation functions involved. After that we 
will discuss the parametrizations, such that we are able to address the Sivers effect, the Collins effect and the 
Qiu-Sterman mechanism (gluonic poles), explicitly. 



A. Distribution functions 



Hermiticity, parity and time reversal invariance yield conditions for the correlator that constrain its 
parametrization, 

$t(p; P, S) = 70 $(p; P, S) 7o [Hermiticity], (91) 

*(p, P, S) = 70 $(p; P, -S) 70 [Parity], (92) 

$*(p; P, S) = (-i 75 C) <P(p; P, S) H 75 C) [Time reversal], (93) 

where C — «7 2 7o, — i 7 5C= i 7 7 and p = (p ,— p). Similar conditions arise for the fragmentation matrix 
elements. Including link operators and for quark-gluon matrix elements slightly different conditions apply. For 
the gauge link one has 

U{ a , S] = VU [a ^ = U [a $, TI7 [0i£] Tt = ^_ Bi _ a , (94) 

for which we used = A^, Pi^^pt = A^) and T A M (£)Tt = This means that the space- 

reversed (time-reversed) correlation function has a different link structure, namely a link running from a (—a) 
respectively. However, if the common point is defined with respect to the two fields in the matrix element, no 



20 



problem arises. For example the straight line link with path z^(s) = (1 — s) M + s£ M gives a path z^ after 
applying parity, but after a change of variables one ends up with the same path; similarly for time reversal. 

For the transverse momentum dependent functions with links running along the minus direction connected 
at infinity, the situation is different. The point n~ = oo is defined by n • n + = 00, which after parity transforms 
into the point fj ■ n + = 00, but after time reversal transforms into the point fj ■ h + = — 00. As a consequence 
one finds that for the p~ -integrated functions 

&+U(x,p T ) = 7o $ [+1 (2^)70 (95) 
<£ l+] (x,p T ) = 7 o$ [+1 (.T,-p T ) 7o (96) 
<&[+]* (x,p T ) = (-i l5 Q&-\x,-p T )(-i l5 C), (97) 

where $H is defined with the link running via n~ = —00, referred to as timelike distribution in Fig. 5. 
Concluding, in the parametrization of & + }(x,p T ) time reversal does not pose constraints. Application of this 
operation transforms <1>H into and vice versa. T-odd quantities will be defined as the ones that vanish 
when $H = $[+1. Accounting for the transformation in Dirac space and the sign change of p T , we define 

2$ [T - cvcn] (x,p T ) = <t> [+] (x,p T ) + (-H 5 C) <P^*(x, -p T ) (-i l5 C) 

= &M(x,p T ) + &-l(x,p T ), (98) 

2& T -° dd \x,p T ) = $M(x,p T ) - {-i l5 C)^*{x,-p T ) (-iTsC) 

= <S> [+] {x iPt ) - <5> [ - ] (x,p T ). (99) 

Note that the name 'T-odd' does not imply a violation of time reversal invariance. For the integrated distri- 
butions the different links ([±]) merge into one, $(x) = & + \x) = $H(:r) and <&(x) is T-even. For transverse 
momentum dependent correlations the sum of $M and <j>H, i.e. <&q(x), is T-even, while the difference, related 
to $q(x,x) (see Eq. 42), is T-odd. Summarizing, we have 

$ [±] {x,p T ) = & T - cvcn \x,p T ) ± & T -° dd \x,p T ), (100) 

and for the integrated and weighted distribution correlators 

$(a;) =$[ T - even l(a;), (101) 
$ [ * ]a (x) =$% [T - even] (x)±iT$% [T - odd] (x,x). (102) 

These results imply that T-odd distribution functions, e.g. the Sivcrs effect appearing in single spin azimuthal 
asymmetries in leptoproduction and in Drell-Yan scattering (^g - '") are opposite in sign [8]. As shown 

in Rcf. [29] the T-odd functions can be considered as imaginary parts in a helicity matrix representation, leading 
to the representation in Fig. 8a. The behavior under time reversal of gluonic matrix elements, like <£>a, ®d and 
<E>G) can also be studied separately It turns out that the quantity $d(x) is T-even, while &g(x, x) is T-odd. 
The latter can also be seen in A + = gauge. Using relation Eq. A8 for £ = — oo~ in the matrix elements 
in Eq. 29 and 30 yields 2n x) — &a(oo)( x > x ) ~ ^A(-oo) ( x > x )> tne gluonic pole matrix element is the 

difference of the boundary terms that transform into each other under time reversal. 



B. Fragmentation functions 

Constraints on the correlator A come from hermiticity, parity and time reversal invariance. The essential 
difference with distribution functions is that time reversal transforms the out-states in the definition of frag- 
mentation matrix elements into in-states. Taking this into account and explicitly adding subscripts in and out, 
one obtains the conditions, 

AH t (z,fc T )= 7o AH(z,/c T ) 7o (103) 
A [ -J t (z, k T ) = 7o A';| (z, -fc T )7o (104) 
A [ -J;(z, k T ) = H75C) Aj+ ] (z, -k T ) H75C), (105) 

where A^ are the spacelike/timelike fragmentation functions, illustrated in Fig. 5. Defining Ao and Apsi as 
sum and differences of matrix elements with out and in-states respectively, one has 
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TABLE II. Summary of time reversal allowed (yes) and forbidden (no) parts in integrated or weighted distribution 
and fragmentation correlators. 





T-even 


T-odd 




T-even 


T-odd 




T-even 


T-odd 




T-even 


T-odd 




yes 


no 


A 


yes 


yes 


Ao 


yes 


no 


Afsi 


no 


yes 




yes 


no 


A e 


yes 


yes 


A SO 


yes 


no 


Aapsi 


no 


yes 


$ D 


yes 


no 


A_D 


yes 


yes 


A DO 


yes 


no 


Adfsi 


no 


yes 




yes 


no 


A A 


yes 


yes 


Aao 


yes 


no 


Aafsi 


no 


yes 


<£>G 


no 


yes 


A g 


yes 


yes 


Ago 


no 


yes 


Agfsi 


yes 


no 



4A^ cvcn, ( 2 , M = AW(z, k T ) + AH(z, k T ) + a{+ ] (z, k T ) + A{; ] (z, k T ), 
4A^- odd] (z,k T ) = A^ t (z,k T ) - A [ -J t (z,k T ) + A^(z,k T ) - A^(z,k T ), 
4Ag- odd] (z,k T ) = A [ +\(z,k T ) + A [ - u \(z,k T ) - A^(z,k T ) A^(z,k T ), 
4 A[, T s 7 cvonl (z, k T ) = AW(z, k T ) - A" (z, k T ) - A [ +\z, k T ) + Aj- ] (z, k T ), 

which implies 

A^(z,k T ) = [A[J- ovcnl (z,fc T ) + A[ T s 7 oddl (z,fc T )] ± [A^- odd] (z,k T ) + Ag- CVCn] ^k T ) 
and for the integrated and weighted fragmentation functions 



A out (z) = A[J- cvcn] (z) 

ol [T — even] 



A£r ddl w, 



\A a 11 



(z) 



A^ T s r ddl (^)] ± [nA a G V-° d %,z) + *A a J T F -r»Hz,z)] . 



(106) 
(107) 
(108) 
(109) 

(110) 



(111) 
(112) 



7T A° 



The essential difference with the distribution functions is that for the fragmentation functions the differences 
between in and out states become relevant. These constitute final state interactions within the soft part (here 
labelled by 'FSF), decoupled from the quark and gluon operators that make the connection to the hard scattering 
part. The effects arising from the difference between [+] and [— ] are labelled by a subscript 'O'. Note that in 
the literature [30,7,15] this is also referred to as initial or final state interactions, depending on the process 
under consideration. The behavior of the various correlators is given in Tabic II. Therefore, in contrast to the 
distribution functions, Ag and Ag contain T-even and T-odd parts. Also the correlators Ajj and A^ contain 
T-even and T-odd parts. 

At first sight Eqs 110 and 112 seem to imply a breaking of universality. Particular transverse momentum 
dependent functions obtained as Dirac projections of A^j.(^, k T ) give unrelated T-odd (and also T-even) results 

for the two possible link configurations, indicated with superscripts [+] or [— ]. In Ag^" t there arise T-odd parts 
from Ag FSI and from 7rA GO , the sign in front of the latter coupled to the link structure. Thus, the T-odd 
(Collins) effects in pion leptoproduction and in electron-positron annihilation are a priori not identical. The fact 
that a certain azimuthal asymmetry arises from a combination of correlation functions (Ag and Ac), however, 
need not imply a breaking of universality. It is quite similar to the q T -integrated structure functions of different 
processes which involve different flavor weights. It seems possible that a factorization proof can be established 
for the correlation functions Ag and Aq separately. 

In Figs 8 we illustrate the differences between and those between the A; n / out . For the former one can 
argue that forgetting about the effect from in- and out-states one has the same situation as for distribution 
functions where the T-odd parts can be considered as imaginary parts of helicity amplitudes (shown as the 
projections along the vertical axis in Fig. 8a). For A in / out one can use the fact that in- and out-states can be 
obtained by different Moller operators, allowing a connection of Aj n and A ou t via a unitary operation, illustrated 
in Fig. 8b. The combined effect is illustrated in Fig. 8c. It shows that the T-odd parts of A^ out 



' + ' in e+e and 



d out 



in SIDIS are in general not equal in magnitude. 



Actually, not only the T-odd effects acquire contributions from both terms in A^ ± Aq, but also the T- 
even effects, hence affecting all comparisons of azimuthal asymmetries in leptoproduction and electron-positron 
annihilation. 
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(a) 



(b) 



(c) 



FIG. 8. Illustration of the role of gauge links ([+] versus [— ]) and the role of final states ([in] versus [out]) in T-odd 
effects for distribution and fragmentation functions. 



Here we will discuss the parametrizations for transverse momentum dependent distribution functions and 
fragmentation functions. The parametrizations of <&(x,p T ) and A(z, k T ) consistent with the conditions imposed 
by hermiticity and parity, including the parts proportional to (M/P + )° and (M/P + ) 1 are given in Appendix B 
[19,6]. In principle the functions could differ depending on the ± gauge link structure for both distribution and 
fragmentation functions and also on the O/FSI characterization for fragmentation functions. 

We will first discuss the parametrizations for the transverse momentum integrated and then the once- weighted 
functions (transverse moments). As explained in Section IV-F, we do not address twice- weighted functions 
(relevant for k T broadening and the average k\ in jets) in view of potential problems with factorization. 

The result for the correlator <!> after integration is 



For the distribution functions the T-odd part vanishes, i.e. ei,{x) = fr(x) — h(x) = 0. For the parametrization 
of fragmentation functions one has both T-even and T-odd functions, i.e. the functions El(z), Dt{z) and H{z) 
appear. Explicitly, 



sub- leading parts of the cross sections. 

Next we consider the p T -weighted results referred to as transverse moments. Restricting ourselves to the 
leading (M/P + )° part, we write for the correlators $g'^ (Eq. 38) 



C. Parametrization of quark and quark-gluon correlation functions 




(113) 
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-tt 1)i±] (x) eV,y»^ - h t m± \x) ^i^ 1 }, (115) 
where we have defined p 2 /2M 2 -moments (transverse moments) as 

Sit 0*0 = J d V 9it(x,P 2 t ), (116) 

and similarly for the other functions. 

For the transverse moments, the average of the ± correlators in Eq. 43 is T-evcn, i.e. 

*g(*) - y { 9 ^ {x) s " 75 ^ + ~ Sl h ^ ){x) 75[7 2 ,A+] }' (117) 

while the gluonic pole contribution &q(x,x) in Eq. 42 is T-odd. Writing 

n $« G (x,x) = ™ I (*) e« - hf W (x) *^±1 1 , (118) 



we have the relations 



5 « [±1 (x)= 3 «(x), (119) 
^i 1)[±1 («) = ^i 1) (*), (120) 



/l ^)W(,) = ± /^)(,), (121) 
ht W[±] (x) = ±ht {1) (x). (122) 



These results show that azimuthal spin asymmetries involving the distribution functions g^ and h^f\ are process 

independent, while those involving the functions and h^ 1 " 1 change sign. Eq. 121 represents the explicit 

connection between the Sivers effect (l.h.s.) and the Qiu-Sterman effect (r.h.s., cf. Ref. [26]); Eq. 122 is its 

chiral-odd counterpart. We note that due to the presence of gluonic pole effects the evolution equations of f^^ 1 
and h^ 1 ^ [41] need to be reconsidered. These functions originate solely from gluonic pole effects. 

For the fragmentation functions, not only the parametrizatio 
parts, but also the average contains T-even and T-odd parts, i.e. 

Ag(z) - M h \zG$(z) S« 75 A- - S hL zH^\z) 



For the fragmentation functions, not only the parametrization for A^' contains both T-evcn and T-odd 



■^W^pf^ - zHt<%)*£M\, (123) 



For the gluonic pole contribution one obtains 

tt Ag(z, z) = M h j*G#(s) S% T 75 fi- - S hL zH±£\z) 75 [7 °'^ ] 

e% vp rn^ hT - ^(,) il^lj. (124) 

and we obtain for the transverse moment in A^^ , 

G^ [±1 (z) = G ( $(z) ± (125) 

^ 1)[±1 (,) = < (1) (,)±^ (1) (,), (126) 

C^ 1)[±1 (z) = £>^(z) ± D±£\z), (127) 

^ (1)[±1 (z) =i/ 1 ±(1) (z)±iJ 1 ±(1) (z). (128) 
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The occurrence of out-states in the fragmentation matrix elements is responsible for the appearance of the 
T-even functions G^ and H^jp and the T-odd functions and H^ X \ We note that these results differ 

from Ref. [19]. For example, Eq. 128 shows the explicit forms of the Collins function in e + e~ (plus sign) 
and SIDIS (minus sign), respectively. Furthermore, the above results imply that the evolution equations of 
G^, H^jp , D^ 1 ^ and H^ 1 * 1 [41] also need to be reconsidered. 

For functions weighted twice with a transverse momentum in &qq, one needs higher transverse moments of 

-L( 2) 

the functions, such as (x). Relations for these functions involve not only twist-three, but also twist-four 
parts in the correlators. This includes the simple pj, average in $,92(2;). As mentioned above, for these functions 
we do not expect a simple process dependence as given for the once-weighted results, but this requires further 
investigation. 

In Rcfs. [42,26] parametrizations have been given for two-argument quark-gluon correlations. As shown in 
Eq. 41 one finds after integration correlation functions with D a tp(x) at the same point, for which the QCD 
equations of motion can be used. These relate <&d{x) (see Eq. 41) and also 70 $d{x) 70 to $(x); explicitly, 

= Y {(*9r(x) + S L (xh L (x) -^9i(xj) 75[7 ; A+] 

-(xe { x)-^ h{ x))^M^ (129) 

which contains only T-even functions. For fragmentation functions one also finds T-odd functions, 

7 5 [7 Q ,A-] 



A£(*) = M h I (g t (z) zH^ S% T 75 ^ + S hL (h l (z) zG 1 (z)^j 



- { E(z)-^-zD l{ z^ 



M h ^ 7 4 

+D T (z) z% vp ^r?_S» hT + H(z) ^llLI + El (z) ^i^zl j. (130) 

Besides the relations resulting from the equations of motion, also Lorentz invariance may lead to relations 
between correlators [43,19,41]. As can be seen from the explicit treatment in Ref. [44], these relations are 
derived from the Lorentz structure of non- integrated quark-quark correlators as in Eq. 9. At present it is not 
clear how matrix elements of A T fields at infinity and hence the link structure play a role in these relations (see 
also Ref. [45]). 



VIII. SUMMARY AND CONCLUSIONS 



In this paper we have analyzed transverse momentum dependent distribution and fragmentation functions 
appearing in several hard processes in which at least two hadrons are involved in initial and/or final state. In 
these processes one has besides a hard scale Q, a non-collinearity q T which is characterized by a hadronic scale 
Q T and an azimuthal angle. We have shown explicitly, using the results of Belitsky et al. [14] how quark-quark- 
gluon matrix elements lead to fully color gauge invariant definitions of the correlation functions that appear 
in leading and first sub- leading order in 1/Q in the hadron tensor of hard processes. The gluon fields appear 
in the gauge link connecting the quark fields. The transverse gluons needed in the gauge link for transverse 
momentum dependent functions involve gluon fields at lightlike infinity. The fact that the gluonic effects can 
be cast into (conjugate) links attached to the two (conjugate) quark fields still allows for an interpretation of 
these functions as probability densities. 

The structure of the gauge links in hard processes is not always the same. In particular the gauge links in 
distribution functions in SIDIS and the DY process run in opposite lightlike directions indicated with indices 
±. The two different correlators are connected via a time reversal operation. Similarly, the gauge links in 
fragmentation functions in SIDIS and electron-positron annihilation run in opposite lightlike directions. At 
leading order, the difference between ± correlators vanishes upon integration over transverse momenta. In q T - 
weighted cross sections, one finds correlators weighted with transverse momentum (transverse moments), which 
are dependent on the (±) link structure. The same quantities appear in subleading integrated cross sections. The 
difference between transverse moments with different (±) link structure corresponds to a (color gauge invariant) 
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gluonic pole matrix element, the word pole referring to the fact that one deals with a zero-momentum gluon 
field. This matrix element appears in different processes as the Qiu-Sterman effect. 

Considering the behavior under time reversal for distribution functions, it turns out that the gluonic pole 
contributions coincide with T-odd effects, leading to single spin asymmetries, like the Sivers effect. This es- 
tablishes the direct connection between the Sivers effect and the Qiu-Sterman effect. Since for distribution 
functions gluonic poles are the only source of T-odd effects, one finds that these effects have opposite signs in 
SIDIS and the DY process. For fragmentation functions, T-odd effects leading to single spin asymmetries, like 
the Collins effect in pion leptoproduction, arise not only from the gluonic pole contribution, but also from final 
state interactions. The latter are purely soft interactions in the fragmentation part and has the same sign in 
different processes. Hence the T-odd effects in SIDIS and electron-positron annihilation are not connected by a 
simple sign relation. This is not a breaking of universality, but rather the appearance of different combinations 
of fragmentation functions in different processes. The T-odd effects will not only appear in Collins asymmetries 
in SIDIS or electron-positron annihilation, but likely also in other processes. We emphasize the importance of 
an analysis of the link structure in processes like pp^ — > -k X . We note that also for distribution functions T-odd 
contributions with the same sign in different processes could arise if time reversal is realized in a nonstandard 
way [46]. This would spoil the simple sign switches in Eqs 121 and 122. 

Gluonic pole contributions appear in the azimuthal asymmetries, in the processes that we have considered, 
with a particular sign. This affects not only T-odd, but (in case of fragmentation) also T-even azimuthal 
asymmetries. The transverse moments with a different link structure differ by 'effective' twist-three functions. 
The consequences of these gluonic pole contributions on the evolution of the transverse moments have not been 
considered. Given the known operator structure of the contributions, however, such a study ought to be doable. 

Although in trying to model distribution or fragmentation functions [44,47] one is always stuck with the 
problem of evolution, modelling has been proven useful to illustrate several effects [48,7,15,49], such as the sign 
change in the Sivers functions and the sign behavior for the Collins function ^\ As seen from Eq. 128, 

the latter can have two contributions with different signs. It is not clear whether both contributions are present 
in the models studied. 

Our general analysis of the various correlators in leading and subleading order single spin and azimuthal 
asymmetries in hard processes may help to analyze the results of the first generation of experiments that 
presently are being performed or planned by, for instance, HERMES (DESY), COMPASS (CERN), BELLE 
(KEK) and the results that are obtained by looking at existing LEP data. We hope to see the emergence of a 
coherent picture of these asymmetries. 
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APPENDIX A: TRANSVERSE GLUON FIELDS AND THE FIELD STRENGTH TENSOR 

In lightcone gauge A + = 0, the relation between A" and G +a becomes 



which can be inverted to yield A" in terms of a boundary term and an integral along the minus direction with 
G +a in the integrand. Without gauge choice we have 



as our starting point. Next we multiply from left and right with link operators U^ a ^ and U^ x a j respectively, 
built from A + fields, running along the minus direction. They are denoted by 
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[ &X ' Arp (X 



')] =G +a (x) 



(Al) 



igG +a {x)= [iD+(x),iD a (x)] = [iD + (x),gA%(x)] - ig [d%, A+(x)] 



(A2) 




(A3) 



and satisfy 
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9iU^ x] =U^ x] D+(x). (A4) 

We then obtain 

i9U^ x] G+»{x)U^ a] = [U^ x] iD+(x)U- ia] ,U- !x] gA^x)U-^ -igU^ [t%,A+(x)] 

= [idi^ x] gA«{x)U^ a] ]-i 9 U^ x] [d% , A+ (x)] U^. (A5) 

Thus we find 

[d + ^ x] A«{x)U^ a] ] = U Wx] (G+»(x)+[dZ,A+(x)}) U { ~ M , (A6) 

which is the relation needed to express the transverse gluon fields in terms of the field strength. In particular 
one has 

^.^(OEfcoo] - ^(°0 = t d T U ^.:n\ ( G+a ^ + l^ A+ (V-)]) t^oop (A7) 

J oo 

or in A + = gauge 

A%(0 - A?(oo") = f dq- G +a ( V ). (A8) 

J OO 



APPENDIX B: PARAMETRIZATIONS OF TRANSVERSE MOMENTUM DEPENDENT 

FUNCTIONS 

The parametrization of &(x,p T ) for a spin or spin 1/2 target, consistent with the conditions imposed by 
hermiticity and parity, including the parts proportional to (M/P + )° and (M/P + ) 1 is given by 

®(x,p T ) = ^-\f 1 (x,pl)fi++g ls (x,p T )'y 5 ifi + 



2 



2 ^ U(x, pi) + f 1 - (x, pl)^+ g ' T (x,pl) 75 p T 



\ 

-Mx,pl) e? lp S Ta S L ft{x,pl) ^-g^ 



Here the spin vector is defined 



e s (x,p T )i 75 + h(x,pl) L '+"- J . (Bl) 



M ^ M 

(5 = for spin 0) and we have used the shorthand notation 



9n(x,p T ) = S L g 1L (x,pl) +g 1T (x,p 2 T ) ^ , (B3) 
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and similarly for h± s , and h s . We note that all noncontracted p T -dependence (including appearance of dot 
products like p T ■ S T ) is treated explicitly, leaving functions depending onpj, which is important in distinguishing 
T-even and T-odd behavior. The structures multiplying the 'twist- two' functions f\, gi s , h\T, h^ s or the 'twist- 
three' functions e, f^, g' T , g^~, h^, h s are T-even (satisfying Eq. 97). The structures multiplying the 'twist- two' 
functions f^ Tl or the 'twist-three' functions fx, , e s , h are T-odd (satisfying Eq. 97 with an additional 
minus sign). 

As notation in the parametrizations for fragmentation we employ capital letters, to be precise 
A(z,k T ) — zDi(z 7 —zk T )fi-+zGi s (z 7 —zk T )j 5 fi- 

+zH 1T (z, ~zk T ) h zH u {z, -zk T ) — — 



+ z£(z, -zfc T ) + zD ± (z, -zk T ) + zG' T (z, -zk T ) 75 fi hT 

P h { M h 

+zGi(z, -zk T ) + zH -(z, -zk T ) + zHs{z , _ zfer) l^M 

M h ZM h Z 

+zD T (z, -zk T ) eP T ° lp S hTa + S hL zD£(z, -zk T ) cg ^ Ta 

Mh 

-zE s {z, -zk T ) J75 + zH{z, -zk T ) lM^i±l|. (B4) 
Here the spin vector is defined 

S h = S h J^-n--?h-n + )+S hT , (B5) 



M h M 

(Sh = for spin 0) and we have used the shorthand notation G\ s , etc., 

(k T ■ Sfir) 

Gis (2, -zk T ) = S hL G 1L (z, -zk T ) + Gi T (z, -zk T ) — T T . (B6) 

Mh 

The second argument of the fragmentation functions is chosen to be k' T = —zk T , which is the transverse 
momentum of the produced hadron with respect to the quark in a frame in which the quark does not have 
transverse momentum. In fact the functions only depend on kr' 2 . As for the distribution functions a division 
can be made into T-even functions (D u G u , H 1T , H^ s , E, , G' T , G^ , H£ , H„) and T-odd functions (D^ T , 
Hj-,D T ,Di,E s ,H). 
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